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In the framework of semihard (kr factorization) QCD approach, we consider the differential cross sections of 
D meson production at HERA. The consideration is based on BFKL and CCFM gluon distributions. We find 
that in the case of BFKL LO gluon distribution the theoretical results are sensitive to the Pomeron intercept 
parameter A. We present a comparison of the theoretical results with available ZEUS experimental data. 



1. INTRODUCTION 



grated out: 



The experimental results on deep inelastic 
charm production obtained by the HI Jl|] and 
ZEUS @,|§ Collaborations at HERA provide a 
strong impetus for further theoretical studies. 
This process is truly semihard because of the pres- 
ence of two larges scales: the virtuality of the ex- 
changed photon (Q 2 ) and the charm mass (m 2 ), 
both being much larger than Aqcd but much 
smaller than s. Therefore, in the present note, 
we focus on the semihard approach |^,|| (SHA), 
which we had applied earlier to the D*^ meson 
photoproduction MM in a similar manner. 



2. THE SEMIHARD APPROACH 



It is known that the resummation y,g|,g|]9j of 

the terms [ln(^ 2 /A 2 ) a s } n , [ln{l/x) a s } n and 
{\n(n 2 /A 2 )ln(l/x)a s ] n in SHA results in the so 
called unintegrated gluon distribution 
F(x, k 2 ^, Qq), which determines the probability to 
find a gluon carrying the longitudinal momen- 
tum fraction x and transverse momentum kx at 
the probing scale Qq. It obeys the BFKL equa- 
tion jl(| and reduces to the conventional gluon 
density G(x,[i 2 ) once the qr dependence is inte- 



3~ (X , krp , 



dk\.=x G(x, ji 2 ). 



(1) 



The factorization scale Qq (such that a s (Ql) <1) 
indicates the scale of the nonperturbative input 
distribution. 

The CCFM evolution equation jllj includes co- 
herence effects via parton angular ordering and 
reproduces the BFKL evolution equation in the 
small x limit. Therefore the CCFM unintegrated 
parton distribution A(x, fc§,, Qq, q 2 ) (unlike the 
function F(x, k\, Qq)) depends also on the max- 
imum angle allowed for any emission correspond- 
ing to q = ptjjj- — z). In the small x limit it 



reduces to T 11 



When calculating the spin average of the ma- 
trix element squared, we substitute the full lepton 
tensor for the photon polarization matrix: 



t 7 t 7 



^ - A(p e q)g^]/(q 



2\2 



(2) 



(including also the photon propagator factor). 
The virtual gluon polarization matrix is taken in 
the form 101: 



p^x 2 /\k T \ 2 = k»k v T /\k T \ 



(3) 



Figure 1. Differential cross sections for deep inelastic D production in the ZEUS accessible kinematical 
region as functions of: (a) log 10 Q 2 , (b) log 10 x, (c) W, (d) pr{D*), (e) rj(D*) and (/) z(D*). 



where p e and p p are the 4-momenta of the incom- 
ing electron and proton. 

To parametrise the unintegrated structure 
functions, we use the prescriptions of ref. (l2| . 
The proposed method lies upon a strai ghtfor- 
ward perturbative solution of the BFKL equa- 
tion where the collinear gluon density xG(x,(jl 2 ) 
is used as the boundary condition in the in- 
tegral form (1). Technically, the unintegrated 
gluon density is calculated as a convolution of the 
collinear gluon density with universal weight fac- 
tors: 

F(x,k 2 t ,n 2 ) = /' 1 g(r ? ,fe t 2 , M 2 )-G(-,M 2 )rfr ? , (4) 
Jx V V 

g( V , kl^) = Jo(2^a s ln(l/77)ln(^/A: 2 )), 
k 2 t < M 2 , (5) 



0fa,fc?,/i 2 ) = ^/ (2V / « S ln(lA7)ln(fc 2 /M 2 )), 

k\ > M 2 , (6) 

where Jo and Iq stand for Bessel functions (of 
real and imaginary arguments, respectively), and 
a s — a s /3n. The latter parameter is connected 
with the Pomeron trajectory intercept: A = 
a s 41n2in the LO, and A = a s A\n2-Na 2 in the 
NLO approximations, respectively, where N ~ 18 

In the previous wor k [p| we used the standard 
GRV parametrization [[L4j for the collinear gluon 
density, from which the unintegrated gluon distri- 
bution was developed according to eqs. (4) - (6). 
Some other essential parameters were chosen as 
follows: the charm quark mass m c = 1.5 GeV, the 
Peterson fragmentation parameter e = 0.06, the 
overall c — > D* fragmentation probability 0.26. 



The Pomeron intercept A was regarded as free 
parameter, and then the value A = 0.35 has been 
extracted from a fit to the experimental pr(D*) 
spectrum measured by the ZEUS collaboration 

In this paper we also use an unintegrated gluon 
density coming from a solution of the CCFM evo- 
lution equation (see [^6|). It was shown in Jl6| 
that a good description of the inclusive structure 
function F%{x, Q 2 ) and the production of forward 
jets in DIS, which arc believed to be a prominent 
signature of small x parton dynamics, can be 
obtained using the CCFM unintegrated gluon 
distribution. In we have used the hadron level 
Monte Carlo program CASCADE described in 
pG[ to predict the cross section for D* photopro- 
duction at HERA energies. It was shown that in 
this case the description of the differential cross 
section da/dr](D*) for Q 2 < 1 GcV 2 for different 
regions of pt(D*) is improved. 
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3. NUMERICAL RESULTS 

The theoretical predictions on the differential 
cross sections of deep inelastic D* production are 
shown in Figure Q for the ZEUS kinematical re- 
gion: 1 < Q 2 < 600 GeV 2 , 1.5 < priD^) < 
15 GeV and r){D*±) < 1.5. Different curves 
in Figure [l] correspond to different values of the 
Pomeron intercept parameter: A = 0.166 (dot- 
ted), 0.35 (solid) and 0.53 (dashed). We see that 
the theoretical curves with A = 0.35 describe 
all ZEUS experimental data || on the differen- 
tial cross sections, except the da/dr](D*) distri- 
bution. The SHA calculations with BFKL un- 
integrated gluon distribution show some shift to 
negative rj{D*) with respect to the data. This 
descrepancy between the data and the SHA pre- 
diction could result from the use of the Peterson 
fragmentation function. Another reason may be 
connected with the BFKL LO unitegrated gluon 
density. 

The use of CCFM unintegrated gluon den- 
sity (with angular ordering) from MC generator 
CASCADE ||| with the Peterson fragmenta- 
tion function for c — > D* + transition gives some 
positive shift to the da/drj(D*) distribution (Fig- 



Figure 2. da/dr](D*), CCFM scheme with Peter- 
son fragmentation. 
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Figure 3. da/dr]{D*), CCFM sheme with JET- 
SET fragmentation. 



ure ||) compared to BFKL curve (see Figure |l]e) . 
Finally, when we use the CCFM unintegrated 
gluon distribution with more realistic JETSET 
based fragmentation function |17| implemented 
in (ilf | we obtain good agreement between our 
theoretical results and the ZEUS experimental 
data § for da/dr](D*) (Figure |). 
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